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LIKELIHOOD-BASED GEOLOCATION PREDICTION ALGORITHMS 
FOR CDMA SYSTEMS USING PILOT STRENGTH MEASUREMENTS 

Cross Reference To Related Applications 

This application is related to U.S. Serial No. 09/139,107 entitled "Pattern 
Recognition-Based Geolocation", filed in the names of T.C. Chiang et al on August 26, 
1998; U.S. Serial No. 09/294,997 entitled "A Bayesian-Update Based Location Prediction 
Method For CDMA Systems", filed in the names of K.K. Chang et al on April 20, 1999; and 

U.S. Serial No. ? (Lucent file #CHANG 4-2-16) entitled "Geolocation Estimation 

Method For CDMA Terminals Based On Pilot Strength Measurements", filed in the names 
of K.K. CHANG et al on May 28, 1999. These related applications are assigned to the 
assignee of the present invention and are meant to be incorporated herein by reference. 

Background Of The Invention 

Field of the Invention 

The present invention relates to a method of locating a mobile telephone unit within 
a cellular service area, and more particularly to a method of predicting the location of a 
CDMA mobile unit based upon the probability of its being at a particular location of the 
service area using an algorithm providing a likelihood estimation of the mobile unit's 
location in response to a sequential set of attributes observed by the mobile unit and reported 
back to a base station. 
Description of Related Art 

A cellular telephone system must be able to locate a mobile unit within a cellular 
service area under various RF propagation conditions such, for example, when an E91 1 call 
is made from the mobile unit. Conventional methods for locating a mobile unit are typically 
based on either a triangulation technique which requires signals from three or more base 
stations within a designated service area, or an angle of arrival technique which requires at 
least two base stations. In many areas, the number of base stations the mobile unit can 
detect is less than two. Furthermore, both the triangulation and angle of arrival techniques 
inherently suffer from inaccuracies and signal fading which result from multi-path 
propagation. 
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In the first above-noted related patent application U.S. Serial No. 09/139, 107 entitled 
"Pattern Recognition-Based Geolocation", RF characteristics pertaining to one or more pilot 
signals radiated from a base station and specific to a particular location within the service 
area are detected by a mobile unit and transmitted back to a base station where they are 
5 matched to a known set of RF characteristics and other information obtained from making 
attribute information measurements at all the grid points (sub-cells) in a cellular service area 
and which are then stored in a database located, for example, in a base station server. 

In the second above-noted related patent application U.S. Serial No. 09/294997 
entitled "A Bayesian-Update Based Location Prediction Method For CDMA systems", the 

10 invention is directed to a method of estimating, by a Bayesian probability algorithm, the 
location of a mobile unit in the service area of a CDMA cellular telephone system using a 
model based approach which, among other things, simplifies the generation of a database 
containing a pilot signal visibility probabilities. This eliminates the need for repeated 
attribute measurements at all of the grid points in the service area. 

15 In the third above-noted related patent application U.S. Serial No. „ 

entitled "Geolocation Estimation Method For CDMA Terminals Based On Pilot Strength 
Measurements", the invention is directed to a method of estimating the location of a mobile 
unit in the service area of a CDMA cellular telephone system also using a model based 
approach, but which now eliminates the need for a stored database containing pilot signal 

2 0 visibility probabilities for all of the grid points or sub-cells in the cellular service area. The 

estimation procedure is based entirely on analytical results involving one or more key 
approximations derived, for example, from an integrated model of the wireless 
communications system, its RF environment, and attribute measurement. 

Summary 

25 The subject invention is directed to predicting the location of a mobile wireless 

communication unit in the service area of a CDMA communications system utilizing two 
likelihood functions that define maximum likelihood estimators of the mobile unit's 
location, based on attribute measurements, such as but not limited to pilot signal strength, 
being made at the location of the mobile unit and reported back to a base station. One of the 

3 0 likelihood functions comprises a frequentist likelihood function and the other comprises a 
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Bayesian-modified likelihood function. The likelihood functions are based on the 
assumption that there is an RF model which provides the probability a mobile unit is able to 
detect one or more attributes associated with an arbitrary base station, given it is located at 
an arbitrary location within the service area. The frequentist likelihood assumes the RF 
model provides exact visibility probabilities. In contrast, the Bayesian-modified likelihood 
assumes the RF model only provides reasonable approximations to the true visibility 
probabilities, and uses the approximations to construct a Bayesian prior distribution for the 
true values. Each of the likelihoods can be used in an iterative fashion to produce a 
maximum likelihood estimator for the location of the mobile unit by determining the 
coordinates within the service area which maximize the respective likelihood function. 
Alternatively, or in addition to, each of the likelihoods can be incorporated into a sequential 
Bayesian procedure which outputs a posterior distribution for the location of the mobile unit. 

Brief Description of the Drawings 

Figure 1 is illustrative of a cellular service area divided into a plurality of cells; 

Figure 2 is illustrative of the cells shown in Figure 1 being further divided into sub- 
cells; 

Figure 3 is illustrative of an embodiment of the subject invention; and 

Figure 4 including Figures 4 A-4E comprise flow charts which are illustrative of the 

preferred methods of determining the locality of a mobile unit in accordance with the subject 

invention. 

Detailed Description of the Invention 
Referring now to the drawings and more particularly to Figure 1, the reference 
numeral 10 denotes a service area for a CDMA cellular telephone system partitioned into 
a plurality of contiguous cells 12j,...,12„ . Figure 1 also depicts a plurality of base 
stations 14 13 ...,14 rt located within the service area 10. Also, the service area 10 includes 
at least one mobile switching center (MSC) 16. Typically each of the base stations 
Mj,...,!^ has a sectorized antenna with a distinct pilot signal channel associated with 

each sector. Three sectored antennas are most common. Each sector of the antenna 
serves a corresponding sector of the associated cell. In Figure 1, all of the base stations 
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14,,...,14 B have three sectors each. The three sectors associated with base stations 14j, 
for example, are denoted by the symbols 15 u , 15 u and 15j 3 , respectively. A mobile unit 
20 is shown in sector 15 2 l . 

Figure 2 is illustrative of the cells 12 1 ,...,12 n being further divided into sub-cells 
18 and which are represented by a grid formed by rectilinear grid lines 20 and 22. The 
reference numbers 1, 2, 3 ... 6 of Figure 2 represent individual sub-cells 18j ,...,18 6 , 
respectively. 

Turning attention now to Figure 3, shown thereat is a diagram broadly illustrative 
of the system architecture for determining the location of a mobile unit 20 within the 
service area 10 in accordance with the subject invention. The MSC 16 operates in 
conjunction with the plurality of base stations 14 15 ...,14 n and connects to the local 
telephone system, not shown. A server 22 including digital computer apparatus 23 and 
memory 24, for storing computation procedures, model parameters and system data, are 
typically located at the site of the MSC 16 for purposes which will now be explained. 

In the invention described in the first referenced related application Serial No. 
09/139,107 entitled "Pattern Recognition-Based Geolocation," each sub-cell 1B X 18 rt 
of the service area 10 is identified by a set of observable characteristics which are 
referred to as attributes. Examples of attributes are pilot signal strengths (Ec/Io), phase- 
offsets, angles of arrival, and pilot round trip delays. The invention of Serial No. 
09/139,107 includes a database which contains attribute information which differentiates 
one sub-cell 18 from another and is generated by making a repeated and exhaustive 
survey which involves taking repeated measurements at all the sub-cells 18j ,...,18 rt 
(Figure 2) of the service area 10. 

During the operation phase, after the database has been set up and the location 
service has been deployed, the mobile unit 20 detects and measures attribute values from 
its actual location in sub r cell 18,. and reports them via a message, e.g., a pilot signal 
strength measurement message (PSMM), to the base station(s) 14,,...,14 w (Figure 3), 
which can be one or more of the base stations with which it is in communication. The 
base station(s) forward their respective reported measurements to the geolocation server 
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22. The digital computer apparatus 23 associated with the server 22 statistically 
compares the measured values with the known attribute values stored in the database 
(memory) 24 of all the sub-cells 18 in the service area 10. The sub-cell 18 f whose 
attribute values as stored in the database provide the best match for the measurements 
reported by the mobile unit 20 is considered to be the best estimate of the mobile unit's 
location. 

In the second referenced related application, Serial No. 09/294,997, entitled "A 
Bayesian-Update Based Location Prediction Method For CDMA systems", a database is 
also used to assist the process of location estimation. However, in contrast to the first 
referenced patent application, i.e. Serial No. 09/139,107, it uses a model based approach 
to generate a database containing pilot visibility probabilities for different sub-cells 18 in 
the service area 10. The model-based approach requires that a limited number of pilot 
strength measurements be carried out along a few representative routes in the service area 
10. These measurements are then used to identify the parameters of the model that 
characterizes the service area and its RF environment. Once these parameters are 
identified, simulations are then carried out to populate the database containing the pilot 
visibility probabilities, which are used in the computation of the location distribution of a 
mobile unit requesting location service. An iterative procedure based on a Bayesian 
probability computation is then used to obtain improved estimates of the mobile unit's 
location in response to multiple sets of attribute measurements being reported by the 
mobile unit 20. The model-based approach eliminates the need to carry out extensive 
measurements required by the first named invention, Serial No. 09/139,107. 

In the third referenced related patent application Serial No. , entitled 

"Geolocation Estimation Method For CDMA Terminals Based On Pilot Strength 
Measurements", the model-based approach embodied in Serial No. 09/294,997, "A 
Bayesian-Update Based Location Prediction Method.... 11 to characterize the RF 
environment is used, as is the iterative procedure for computing the Bayesian posterior 
distribution for the location of the mobile. However, the database containing pilot 
visibility probabilities is replaced by analytical formulas that can be evaluated in real 
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time. The evaluation procedures are compact and can typically be evaluated in the digital 
computer apparatus 23 shown in Figure 3. 

Considering the present invention, the analytic formulation for the pilot visibility 

probabilities taught in the third referenced application Serial No. , 

5 "Geolocation Estimation Method For CDMA terminals Based On Pilot Strength 
Measurements", now serve as the starting point for the derivation of two likelihood 
functions, hereafter referred to as the frequentist and Bayes-modified likelihood functions, 
respectively. Each of the likelihood functions is derived based on the assumptions and 
mathematical formulations described in attached Appendix A. In as much as the likelihood 
^ 10 functions depend on the analytic evaluation of the pilot visibility probabilities, attached 
J3 Appendix B provides a self-contained development of the relevant details of these formulas. 

Ill Each likelihood function is a function of (x, y), an arbitrary location of the mobile unit 20. 

^ Accordingly, each likelihood function could be used to obtain a maximum likelihood (ML) 

£ estimator of the location of the mobile unit 20 by finding the (x, y) which maximizes the 

y ~ 15 value of the respective likelihood function. An iterative technique for sequentially updating 
Q each ML estimator with additional pilot signal strength measurements can be utilized. Each 

pi of the two likelihood functions can also be incorporated into a sequential Bayesian 

6j procedure, which outputs a posterior distribution for the location of the mobile, 

yy The Bayes-modified likelihood function, whether it is used in the context of ML 

2 0 estimation or a sequential Bayesian procedure, is a substantial deviation from the inventions 

disclosed in the second and third above referenced related applications in the following way. 
Both of these previously disclosed inventions use an RF model to estimate pilot visibility 
probabilities, the former via simulation techniques, the latter via analytical formula 
evaluation, and implicitly assume that the model holds precisely. The present invention, 
25 however, uses the same RF model only to determine the means of beta distributions that are 
used as Bayesian priors for the true (unknown) pilot visibility probabilities. A beta 
distribution is completely determined once its mean and variance have been specified. 
Subject to the fixed mean values, each of the beta distributions is fully specified by 
maximizing their variances. Maximizing the variance of the beta priors, subject to the 

3 0 specified mean values, is consistent with a non-informative (vague) prior specification. 
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A summary of the derivations which appear in Appendix A will now be given as 
prefatory remarks to the description of the overall mobile unit 20 location prediction process 
depicted in Figure 4. In the example shown in Figures 1 and 2, if the primary base station 
for the mobile unit 20 is 142, the region designated as 17 in Figures 1 and 2 is the set of 
5 feasible locations of the mobile unit 20, and will hereafter be referred to as the region A. 
Along with the region A, the set of all pilots K which are likely to be visible at some of the 
grid points 18 in the set A is defined. 

For each (x,y) <= A , let &#(x 9 y) denote the true probability that the mobile unit 20 

is able to see the pilot in sector j of base station / when it is located at (x y y) . From hereon, 

10 the notation ij will be used to exclusively reference pilots from the set K. Let 0y(x 9 y) 

denote an approximation of 0 v (x 9 y) based on an RF model described in Appendix B. For 

each pilot ij y let jf iS equal one or zero depending on whether the mobile unit 20 can see pilot 

ij at the s-th measurement epoch or not, respectively. The frequentist likelihood through the 
first s measurement epochs has the following recursive form, starting with the definition 
15 L^ix^y)^!: 

£lft(*,jO<*£t^ > (*,y)*A (1) 

For the Bayes-modified likelihood, the prior for 0^ (x,y) is a beta distribution with 
parameters: 



20 



1 Jf 9 i} {x,y)<H2 

- J ^ L - ,i/0 l3 (x,y)>l/2 (2 > 



l-0Jx,y) 



(3) 

1 ,ifd u {x,y)>\/2. 
25 The Bayes-modified likelihood function through the first s measurement epochs has 

the following recursive form starting with the definition I? BlulL {x,y)=\ for all 
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(x,y)eA 



s-\ 

where nf 1 = ^ //* is the number of times pilot ij was visible amongst the first 5-1 

jt=i 

measurement epochs. 

Each of the likelihood functions (1) and (4) are functions of (x 9 y) e A , an arbitrary 
possible location for the mobile unit 20. The ML estimator for the location of the mobile 
unit 20 is obtained by evaluating (1) for all (x 9 y)eA and selecting the values, say 

( x Ll > yldL) > which gives the largest value of (1). An updated ML estimate is produced at 
each measurement epoch. In a similar way, using function (4) rather than function (1) 
generates a sequence of Bayes-modified ML estimates, say (x^ , y%^) . 

Utilizing functions (1) or (4) with a Bayesian sequential procedure can generate an 
alternative sequence of predictions for the location of the mobile unit 20. In each case, the 
initial prior distribution for the location of the mobile unit 20 is assumed to be a discrete 
uniform distribution of the form: 

^(^) = Cl(^) = pj , (x 9 y)sA (5) 

where ||^4|| is the number of grid points 18j,.. AS n contained within A, The posterior 

distribution for the location of the mobile unit 20, through s measurement epochs, based on 
the frequentist likelihood function (1) is, up to a constant of proportionality: 

Alternatively, the posterior distribution of the location of the mobile unit 20, through s 
measurement epochs, based on the Bayes-modified likelihood function (4) is, up to a 
constant of proportionality: 



1 BML 



£J a v (x,y)+0 v (x,y) + s-l 
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A Bayesian sequence of predictions on where the mobile unit 20 is located follows from 
functions (6) or (7) by using the mean or mode of the posterior distribution obtained at each 
measurement epoch. When function (7) is used, the sequence of prediction involves two 
distinct prior distributions, beta and discrete uniform, and the methodology is referred to as 
5 doubly-Bayesian. This completes the summary of Appendix A. 

A description of the mobile unit 20 location prediction process depicted in Figure 4 
will now be given. Considered in light of the accompanying appendices A and B and 
referring to Figures 4A-4E, the location prediction process in accordance with the subject 
invention, as noted above, is implemented in software which resides in the computer 
p% 10 apparatus 23 located at the geolocation server 22 (Figure 3). 

jD The process, referred to hereinafter as the geolocation process, begins at step 30 

yi (Figure 4A) where it is waiting for a new location request to arrive. The geolocation process 

~f continually checks for the arrival of a location request at step 32, and if no such request has 

ap arrived, it goes back to the waiting state (step 30). When a location service request arrives, 

K = S 

15 the geolocation process identifies the domain of support, i.e., the set of feasible locations, for 
the mobile's location based on either the identity of the primary base station of the mobile 
ry unit or the identity of the strongest pilot signal reported by the mobile unit. This is indicated 

%: in step 34. Pilot signals will hereinafter be referred to simply as "pilots." 

yJ In the example shown in Figures 1 and 2, if the primary base station for the mobile 

20 unit is base station^, the region designated 17 in Figures 1 and 2 would be selected as the 
domain of support for the mobile's location and will be hereinafter referred to as the region 
A. Along with the region^, the geolocation process at step 34 also identifies the set of all 
pilots K which are likely to be visible at some of the gridpoints 18 in the region^. Next, in 
step 36 approximations for the conditional probability, conditioned on the mobile unit being 

25 at a location (x,y) gA 9 that each pilot in K is visible to the mobile are computed using 
equation (B9) from Appendix B. At this point, two decisions external to the geolocation 
process must be made. The first decision (step 38) that must be made is which of the ML 
estimation or sequential Bayes estimation methods should be used. For each estimation 
procedure (ML or sequential Bayes) the second decision to be made is which type of 

30 likelihood, frequentist or Bayes-modified, should be used (steps 40 and 42). These two 
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decisions result in the four paths marked (2)-(5). Figures 4B-4E correspond to each of the 
four paths, only one of which would typically be used in any implementation of the 
geolocation process. 

First, suppose path (2) of Figure 4 A is chosen which represents the combination of 
5 the sequential Bayes estimation method and the frequentist likelihood function. This is 
shown in Figure 4B. There step 44 assigns a discrete uniform prior probability to all grid 
points 18 in the set A. The discrete uniform prior reflects the initial state of no information 
about the mobile unit's whereabouts, other than the fact that it resides in the region A. The 
frequentist likelihood function is initialized to unity and time is set to 5=1. Step 46 then 

10 evaluates the frequentist likelihood based on the first set (s=l) of visibility measurements 
using expression (1) noted above. In step 48, the posterior distribution based on the first set 
of visibility measurements is evaluated using function (6). Expression (6) results in values 
that must be normalized so that when they are summed over all (x,y) e A , the result will be 
unity. The posterior distribution gives the updated probability distribution for the location of 

15 the mobile unit 20. A prediction of the location for the mobile unit 20 is next obtained in 
step 50 by computing either the mean or mode of the posterior distribution obtained from 
step 48. The prediction obtained in step 50 corresponds to the first set of measurements 
(s=l). If no further measurements are expected, the geolocation process terminates, 
otherwise as shown in step 52 it proceeds to a waiting state (step 54) and stays there (via 

2 0 step 56) until another set of measurements is received. At that point, the geolocation process 
proceeds to step 58 and increments time up to 5+1 before proceeding back to step 46 and 
looping once again through step 48 and step 50 which lead to an updated prediction for the 
location of the mobile unit 20. Eventually, no further measurements will be expected and 
the geolocation process will terminate at step 52. 

25 Next, suppose path (3) of Figure 4 A, which is depicted in Figure 4C, is chosen 

which represents the combination of the sequential Bayes estimation method and the Bayes- 
modified likelihood function. Step 60 directs the calculation, via equations (2) and (3), 
noted above, of the two parameters for the beta prior distribution that is used for the true 
unknown pilot visibility probabilities. Step 62 assigns a discrete uniform prior probability to 

30 all grid points 18 in the set A. The discrete uniform prior reflects the initial state of no 
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information about the mobile unit's whereabouts, other than the fact that it resides in the 
region A. The Bayes-modified likelihood function is initialized to unity and time is set to 
5=1. Step 64 then evaluates the Bayes-modified likelihood based on the first set (5=1) of 
visibility measurements using functional expression (4). In step 66, the posterior 
5 distribution based on the first set of visibility measurements is evaluated using equation (7). 
Equation (7) gives values that must be normalized so that when they are summed over all 
(x y y) , the result will be unity. The posterior distribution gives the updated probability 
distribution for the location of the mobile unit 20. A prediction of the location for the 
mobile unit 20 is obtained in step 68 by computing either the mean or mode of the posterior 

10 distribution obtained from step 66. The prediction obtained in step 68 corresponds to the 
first set of measurements (5=1). If no further measurements are expected, the geolocation 
process terminates, otherwise as shown in step 70 it proceeds to a waiting state (step 72) and 
stays there (via step 74) until another set of measurements is received. At that point, the 
geolocation process proceeds to step 76 and increments time up to 5+1 before proceeding 

15 back to step 64 and looping once again through step 66 and step 68 which lead to an updated 
prediction for the location of the mobile unit 20. Eventually, no further measurements will 
be expected and the geolocation process will terminate at step 70. 

Next, suppose path (4) of Figure 4 A is chosen which represents the combination of 
the ML estimation method and the Bayes-modified likelihood function. Path (4) is shown in 

20 Figure 4D where step 78 first directs the calculation, via equations (2) and (3), of the two 
parameters for the beta prior distributions that are used for the true unknown pilot visibility 
probabilities. Step 80 initializes the Bayes-modified likelihood function to unity and time is 
set to 5=1. Step 82 then evaluates, for all (x,y) € A , the Bayes-modified likelihood based 
on the first set (5=1) of visibility measurements using expression (4). In step 84, the value of 

25 (x y y)eA which gives the largest likelihood value is selected as the Bayes-modified 
maximum likelihood estimator of the mobile unit's location The estimate obtained in step 
84 corresponds to the first set of measurements (5=1). If no further measurements are 
expected, the geolocation process terminates, otherwise as shown in step 86 it proceeds to a 
waiting state (step 88) and stays there (via step 90) until another set of measurements is 

30 received. At that point, the geolocation process proceeds to step 92 and increments time up 
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to 5^1 before proceeding back to step 82 and step 84 which lead to an updated prediction for 
the location of the mobile unit 20. Eventually, no further measurements will be expected 
and the geolocation process will terminate at step 86. 

Finally, suppose path (5) of Figure 4A and further shown in Figure 4E, is chosen 
5 which represents the combination of the ML estimation method and the frequentist 
likelihood function. As shown, step 94 initializes the frequentist likelihood function to unity 
and time is set to s=l. Step 96 then evaluates, for all (x,y) e A , the frequentist likelihood 
based on the first set (s=l) of visibility measurements using functional expression (1), In 
step 98, the value of (x,y) which gives the largest likelihood value is selected as the 
q 10 maximum likelihood estimator of the mobile unit's location. The estimate obtained in step 
fi 98 corresponds to the first set of measurements (5=1). If no further measurements are 

10 expected, the geolocation process terminates, otherwise as shown in step 100 it proceeds to a 

p waiting state (step 102) and stays there (via step 104) until another set of measurements is 

j; received. At that point, the geolocation process proceeds to step 106 and increments time up 

s 15 to s*l before proceeding back to step 96 and step 98 which lead to an updated prediction for 
the location of the mobile unit 20. Eventually, no further measurements will be expected 
HJ and the geolocation process will terminate at step 100. 

tp Each of the iterative procedures described above with respect to paths (2)-(5), 

^ operates to provide an improved estimate of the location of the mobile unit 20 as more and 

20 more sets of pilot measurements are reported. Each of the procedures can be readily 
extended to include other measured quantities such as phase offsets by including those data 
in the respective likelihood formulations. Also, when desirable, the method of the present 
invention can be modified so that the actual signal strength of the visible pilots can be used 
rather than a binary representation of the fact that the pilots are or are not visible. 

25 The foregoing description of the preferred embodiment has been presented to 

illustrate the invention without intent to be exhaustive or to limit the invention to the form 
disclosed. In applying the invention, modification and variations can made by those skilled 
in the pertaining art without departing from the scope and spirit of the invention. It is 
intended that the scope of the invention be defined by the claims appended hereto, and their 

30 equivalents. 
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Appendix A 

Derivation of the Frequentist and Bayes-modified Likelihood Functions 

5 The notation used in this Appendix is the same as what has been introduced in the 

main portion of this application. In what follows, we make the following assumptions: 

(i) The visibility of a particular pilot at (x,y) e A at a given time t is independent of 
whether the pilot was visible at any time s<t, and is independent of / as well. 

(ii) The visibility of pilot / at time t is independent of the visibility of pilot j (j ^ i) at 

10 time/. 

(iii) Throughout the duration of the prediction process, the mobile is relatively 
stationary, i.e., the mobile unit remains in its initial sub-cell. 

Defining I?(x y y) = 1 , V (x,y) gA , a recursive form for the exact likelihood 
function through the first s measurement epochs for the unknown location of the mobile 
1 5 based on (i)-(iii) is 

L'(x,y)ccL~\x 9 y)^ , (x 9 y)eA. (Al) 

The functions 0^{x 9 y) in (Al) are unknown, so the exact likelihood function is not 
calculable. Replacing the 0 9 (x 9 y) with their approximations 0 ij (x 9 y), defined by equation 

(B9) in Appendix B, gives the definition of the frequentist likelihood function U^fay) 
2 0 shown in expression (1). 

The frequentist likelihood function (1) was derived under the assumption that a 
family of functions 0 jf (x 9 y) accurately approximate the unknown functions 0y(x,y). The 

approximations 0 9 (x 9 y) are based on the RF model described in Appendix B. The RF 
model is an approximation to a very complex stochastic process. Bayesian methods allow 
25 the uncertainty associated with 0 v (x,y) to be incorporated into the model framework. 

Rather than assuming 0 iJ (x,y)zz0 iJ (x,y) , we can model the parameters # iy (x Si y) with 
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independent prior distributions that have mean values equal to 9 tJ {x,y). For each 
(x,y)eA, we use a prior distribution for 9 i} ix,y) that is a beta distribution with 
probability density function given by 

/„, = ^^(i-^^- 1 , 0<u<\ (A2) 

5 where 5[ v ] is the complete beta function. Note that the two parameters, a^x.y) and 

fin ( x >y), to be determined in what follows, depend on both the possible location (x,y) e A 

and the pilot ij sK. Whereas the frequentist likelihood is defined by implicitly assuming 

O Pr[^ iy (x 5 >') = ^(x, i y)]=l, the Bayes-modified likelihood is defined by taking the 

LJ expected value of the exact likelihood (Al) with respect to the beta distributions defined by 

"2 10 (A2). The recursive form of the exact likelihood (Al) does not amend itself to evaluating 
IP the expected value as easily as the non-recursive form of the exact likelihood which is 

k ^(x 5 ^)a:n[^(^^[l"^(^^)r^ , (x>y)^A (A3) 

£ . ye*: 

SJ where is the number of times pilot ij was visible through the first s measurement epochs, 

gr- Evaluating the expected value of (A3), with respect to the beta distributions defined by (A2) 

'% 1 5 gives the Bayes-modified likelihood function 



20 



It is easily verified that equation (A4) has the recursive form shown in the functional 
expression (4). 

25 The definition of the Bayes-modified likelihood function will be complete once 

values for a^x.y) and fi i} (x y y) are specified. Determining what values to use for 



15 

!DS#1 1671 1 /CHANG 5-3-17-19 
Attorney Docket No. 2925-253P 

a ij( x >y) an d Pij(x*y) is an example of the classic Bayesian dilemma - how to specify a 

prior distribution? An underlying principle in specifying a prior distribution is that any 
available information about the parameter of interest, in this case O^^x.y), should be 

reflected in the prior distribution that is used. In this case, we have an approximation of 
5 & iJ (x 9 y) from the RF model, namely 0 lJ (x 9 y) 9 which we will use as the mean of the prior 
distribution for #, y (x,j>) . Accordingly, we have the following constraint: 

Since we have two unknowns, &q(x 9 y) and fi if (x 9 y) 9 we need one additional constraint to 

uniquely determine both a ij {x 9 y) and fi v (x 9 y). To obtain the second constraint, we 

10 impose that the prior be as vague as possible, subject to the first constraint. To implement 
this second constraint, we will maximize the variance of the prior, subject to the fixed mean 
value. The second constraint is consistent with the quantity of the prior knowledge we have 

about 0g(x 9 y) and except for anticipating a value near O v (x 9 y) 9 we know nothing else 

about its potential value. We also restrict attention to beta distributions which have 
15 bounded density functions, implying a iJ {x 9 y)>\ and fi^(x 9 y)^l. The restriction to 

bounded density functions only rules out "U-shaped and "Slide-shaped" beta distributions 
which would not typically be intuitive representations of the prior information on 0 jJ (x 9 y) . 

The problem of selecting a (x 9 y) and J3 tj (x y y) has now been sufficiently constrained so 

that it can be solved by using the following Lemma. 
20 Lemma 

Consider the family of beta distribution with parameters a > 1 and b > 1 . The values of a 
and b which maximize the variance, subject to the distribution having a given mean value 
//,are: a = l and b = (1 -//)///, if //< 1/2; a = ju/(l-ju) and b = 1 , if ju> 1/2. 
Proof of Lemma 

2 5 The variance of the beta distribution is given by a 2 =ab l[(a + b + \){a + b) 2 ]. Since we 
have ju = a/(a+b) 9 wecan write b = a(l -ju)//u 9 and hence a 2 = // 2 (1 - //) / (a + //). To 
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maximize the variance, we need to minimize a, subject to the constraint that both a > 1 and 
b > 1 . Clearly, if ju < 1 / 2 then we can choose a = 1 and b = ju)/ ju>l will be 
ensured. Likewise, if // > 1 / 2 then the smallest we can choose a and still have b > 1 is 
<z = ///(l-//)>l. 

Applying the Lemma with a = a ff (x,^), 6 = ftij( x ,y) and // = 0 ff (x yi y) gives the 
formulas in equations (2) and (3) for a i} (x,y) and fi tJ (x,y) . 
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Appendix B 

Derivation of the Approximations to the Pilot Visibility Probabilities 

In this Appendix, we drive the approximations 0 ff (x,j/) of the pilot visibility 

probabilities. The notation used in this Appendix is the same as what has been introduced 
earlier in this specification. Consider a grid point (x,y) e A . As in the above referenced 

related application Serial No. , entitled, "Geolocation Estimation Method For 

CDMA Terminals Based On Pilot Strength Measurements", the RF power received by the 
mobile unit 20 from sector j of base station / is modeled by the expression 

Rytey) = 2JG,(x, y)I v (x 9 y)F v (x 9 y)M ff (x 9 y) (Bl) 

where 2£ is the transmit power associated with the sector, Gy{x>y) is the antenna gain for 
the sector along the direction pointing towards the location (x 9 y) gA, L tJ (x,y) is the 
distance loss between the base station associated with the sector and the location (x,y) e A , 
Fij( x *y) * s the shadow fading factor and M^(x 9 y) is the measurement noise factor, all in 

absolute, not dB, units. The measurement noise factor is meant to include the effects of fast 
fading (e.g., Rayleigh/Rician) as well as inaccuracies in the measurement process. If 
Y denotes the fraction of T n that is used for the pilot channel, then }R^(x 9 y) is the pilot 

channel power received by the mobile unit 20 when it is located at (x 9 y) e A . 
The model assumes that the distance loss L^{x y y) can be expressed as 

L iJ (x,y) = C p [d ij (x,y)r <B2) 
where d v (x 9 y) is the distance between the base station associated with the sector and the 
location (x,y) e A , and C p and a are constants. Typically, C p takes a value in the range 
10~ 15 to 10" 10 and a is between 3 and 5 when d iJ (x,y)is expressed in miles. The two 
parameters, C p and a , moreover are environment specific. 

The shadow fading factor, F^(x 9 y), models the impact of terrain and large 
structures (e.g., buildings) on signal propagation which create deviations around the signal 



18 

IDS#1 1671 1 /CHANG 5-3-17-19 
Attorney Docket No. 2925-253P 

attenuation predicted by the deterministic path loss factor. We assume that F^{x^y) is 
lognormally distributed and thus can be written as 

F,{x 9 y) = e^e«™ (B3) 
where 0 iy (x, y) is a zero mean Gaussian random variable with standard deviation cr^ . 
5 The measurement noise factor, M^x^y) is also assumed to have a lognormal 

distribution so that we can write 

M^x 9 y) = e&e*> ix ^ (B4) 

q where t$$(x 9 y) is a zero mean Gaussian random variable with standard deviation a M . 

2 Assume now that the pilot strength measurement carried out and reported to the base 

10 station by the mobile unit 20 is the Ec/Io value of the corresponding pilot channel signal. 
01 This value is the ratio of the pilot channel power from the concerned sector received by the 

mobile unit 20 to the total power received by the mobile unit including thermal noise, and 
possibly external interference. 

By letting P^(x 9 y) denote the strength of the pilot channel associated with the 

15 sector as measured by the mobile unit 20 located at (x,y) € A , it follows that 

N 0 + Z R *(*>y) 



Inn 



P v (x,y)= „ (B5) 



where denotes thermal noise and external interference and (as before) kl is shorthand 
notation for the pilot associated with sector / of base station k. 
For convenience, we define 
20 y) = Tip^x^L^y). (B6) 

Observe that the expected value of R^(x 9 y) is equal toQ(x 5t y) . It is implicitly assumed 

that the shadow fading and measurement noise factors are uncorrelated. We approximate 
Py (x 9 y) (B5) by the following expression 
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z ( , Jfy(*>j0 

kleK 
kl*ij 



(B7) 



The difference between (B7) and (B5) is that for kl * ij, the R kl (x,y) quantities have been 
replaced by their mean values. Using (B3) and (B4), we can write 

be* 



a + ce* 



(B8) 



where 



kleK 
kteij 

c = C ij (x,y) 

The probability that the mobile unit receives pilot ij with a signal strength in excess of Tis 

0 ij (x,y) = M^(x,y)>T] 
= Pr[Z l} (x,y)>T] 

- o 9 {x,y) • 

It can be easily shown that 



V t (x,y) = Q 



Ln\^T- 1 -* 

.a a 



(B9) 



where a = (or* + cr* ) 1/2 and Q(u) = -=L [* e~*'* dt . 

V2^ Ju 



Note that the expression for 0 9 (x 9 y) depends on the following unknown 
parameters: a 2 M , N 09 C p and a . As in the above noted related application entitled, 
"Geolocation Estimation in the present invention, pilot strength measurements are first 
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carried out along a small number of representative routes in the service area. The 
measurements so obtained are then used to select values of cr^, cr* , N 0 , C p and a which 

provide the best fit between the predicted pilot visibility probabilities, using equation (B9) 5 
and the observed pilot visibility probabilities 
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Claims 



1 1 . A method of predicting the location of a CDMA mobile communications unit in a 

2 wireless communications service area, comprising the steps of: 

3 (a) receiving measured attribute information from a mobile unit location in the 

4 service area, said attribute information being specific to the location of the mobile unit in the^ 

5 service area; 

6 (b) computing the probability of the mobile unit being at a specific location in the 

7 service area in response to said received attribute information using a likelihood probability 

8 function; 

9 (c) generating an output indicative of the likelihood of the mobile unit being at said 
1 0 location in the service area. 

1 2. A method according to claim 1 wherein said likelihood probability function 

2 comprises an iterative procedure for producing a maximum likelihood estimator of the 

3 mobile unit's location in this service area. 

1 3. A method according to claim 2 wherein the maximum likelihood estimator 

2 comprises the coordinates in the service areawhich maximizes the likelihood probability 

3 function. 



1 4. A method according to claim 2 wherein said likelihood probability function 

2 comprises a frequentist likelihood function. ~" s 

1 5. A method according to claim 2 wherein said likelihood probability function 

2 comprises a Bayes-modified likelihood functtoio 



1 
2 



6. A method according to claim 1 wherein said procedure comprises a sequential 
Bayesian type of procedure. 
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1 7. A method according to claim 6 wherein said sequential Bayesian type of 

2 procedure uses a frequentist likelihood function. 

1 8, A method according to claim 6 wherein said sequential Bayesian type of 

2 procedure uses a Bayes-modified likelihood Sanction. 

1 9. A method according to claim 1 and wherein said attribute information comprises 

2 a pilot signal strength measurement of at least one visible pilot signal at said location of the 

3 mobile unit. ' 

1 10. A method according to claim 9 and additionally including a step (d) prior to step 

2 (b) of reducing the effective size of the service area so as to reduce the number of 

3 computations required in step (b). 

1 11. A method of claim 9 and following step (a) and before step (b), additionally 

2 including a step (d) of identifying a region A of support for the mobile unit in the service 

3 area and (e) identifying a set of possible pilot signals which can be detected by the mobile 

4 unit in said region of support. 

1 12. A method of claim 1 1 and additionally including a step (f) of computing an 

2 approximation 0 fy (x,j>) of the probability the mobile unit detects a pilot signal at said 

3 specific location (x, y) in said region A for all pilot signals in said set of possible pilot 

4 signals and where j is a sector of a multi-sector base station i in said region of support. 

1 13. A method of claim 12 wherein 0^(x 9 y) is derived from an RF model for the RF 

2 power R~ (jc, y) received by the mobile unit and where the model is of the form of 



23 

!DS#1 1671 1 /CHANG 5-3-17-19 
Attorney Docket No. 2925-253P 



3 where T fj is the transmit power associated with the sector j of base station i, (jc, y) is the 

4 antenna gain for the sector j of base station i along the direction pointing towards the 

5 location (x, y) within A, L if (x, y) is the distance loss between the base station i associated 

6 with the sector j and the location (x, y) within A, F (x,y) is the shadow fading factor and 

7 M tJ (x,y) is the measurement noise factor. 

1 14. A method according to claim 13 wherein said attribute information comprises 

2 the measured visibility of a pilot signSltrahsmitted from a base station and the likelihood 

3 probability function comprises an iterative freuentist maximum likelihood (ML) function in 

4 the form of 
5 

1 

8 

9 where x and y comprise rectilinear coordinates in the service area A, s is the number 

10 of measurement epochs, and //£is equal to one or zero depending on whether the mobile 

11 unit can detect a pilot signal ij at measurement epoch s, 

1 15. A method according to claim 13 wherein said attribute information comprises 

2 the measured visibility of a pilot signal transmitted from a base station and the likelihood 

3 probability function comprises an iterative Bayes-modified maximum likelihood (ML) 

4 estimator in the form of 
5 
6 
7 

8 where x and y comprise rectilinear coordinates in the service area A s is the number of 

9 measurement epochs, nff 1 is the number of times a pilot signal in section j of a multi-sector 

10 base station i is visible through the first s-1 measurement epochs, a iJ {x,y) and p i3 {x y y) are 



t*k a 9 (x,y) + fi v (x,y) + s-l 
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11 parameters of a beta distribution, and ju* is equal to one or zero depending on whether the 

12 mobile unit can detect a pilot signal ij at measurement epoch s. 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 



16. A method according to claim 15 where 



1 ( 



,if0 9 (x,y)<l/2 
,i/0 iJ (x,y)>V2 



1 Jf0 iJ (x,y)>l/2. 



1 
2 

3 
4 



6 
7 
8 
9 
10 
11 
12 
13 



17. A method according to claim 13 wherein the frequentist likelihood function 
L^ix^y) is combined with a discrete uniform prior distribution for the location of the 
mobile unit of the form 



14 

to generate a sequential Bayesian procedure which provides a posterior distribution for the 
location of the mobile unit of the form 



where |^4|| is the number of grid points contained within A, and where x and y comprise 
rectilinear coordinates in the service area A, s is the number of measurement epochs, and 
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14 Mij ls equal to one or zero depending on whether the mobile unit can detect a pilot signal 

15 ij.at measurement epoch s. 

1 18. A method according to claim 13 wherein the frequentist likelihood function 

2 L^ix.y) is combined with a discrete uniform prior distribution for the location of the 

3 mobile unit of the form _ 



8 



5 
6 



4 



6 to generate a sequential Bayesian procedure which provides a posterior distribution for the 

7 location of the mobile unit of the form 



11 

12 where \A\ is the number of grid points contained within A, and where x and y comprise 

13 rectilinear coordinates in the service area A, s is the number of measurement epochs, n is the 

14 number of times a pilot signal in section j of a multi-sector base station i is visible through 

15 the first s-1 measurement epochs, a v (x,y) and ft v (x,y) are parameters of a beta 

16 distribution, and fi\ is equal to one or zero depending on whether the mobile unit can detect 

17 a pilot signal ij at measurement epoch s. 

1 19. A method according to claim 1 8 where 

2 



3 a 0 (x,y) = < 

4 



P<M,y) = 



. 1 ,if0„(x,y)£l/2 
6Jx t y) 

, ~/\ Jf0,M,y)>\l2 
l-0 o (x,y) 

\-e y (x,y) 



1 ,if0 y (x,y)> 1/2. 
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1 20. Apparatus for estimating the location of a mobile communications unit in the 

2 service area of a wireless communications system comprising: 

3 a plurality of base stations and at least one switching center providing a common 

4 access to said plurality of base stations; 

5 computing apparatus including a memory, electrically connected to at least one of 

6 said base stations or said switching center, for storing a set of parameters of a model of the 

7 wireless communications system including the RF environment and predicted values of 

8 attribute measurement of the service area and wherein the values of the parameters are 

9 adjusted to obtain a substantial match between measured attribute values at a predetermined 

10 number of locations in the service area and the corresponding initial values predicted by the 

11 model; 

12 said computing apparatus further including software for calculating, in response to 

13 one or more attribute values being measured and reported by the mobile unit from a specific 

14 location within the service area, a predicted location of the mobile unit within the service 

15 area using a likelihood probability function; and, 

16 circuit means coupled to said computing apparatus for generating an output 

17 indicative of the predicted location of the mobile unit. 
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ABSTRACT 



The location of a mobile wireless communication unit in the service area of a 
CDMA communications system is predicted utilizing two likelihood functions that define 
maximum likelihood estimators of the mobile unit's location, based on attribute 
measurements, such as but not limited to pilot signal strength, being made at the location of 
the mobile unit and reported back to a base station. One of the likelihood functions 
comprises a frequentist likelihood function and the other comprises a Bayesian-modified 
likelihood function. The likelihood functions are based on the assumption that there is an 
RF model which provides the probability a mobile unit is able to detect one or more 
attributes associated with an arbitrary base station, given it is located at an arbitrary location 
within the service area. Each of the likelihoods are also incorporated into a sequential 
Bayesian procedure which outputs a posterior distribution indicative of the location of the 
mobile unit. 
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IN THE UNITED STATES 
PATENT AND TRADEMARK OFFICE 

Declaration and Power of Attorney 
As the below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated below next to my name. 

I/We believe I/We am/are the original, inventor(s) of the subject matter which is 
claimed and for which a patent is sought on the invention entitled LIKELIHOOD-BASED 
GEOLOCATION PREDICTION ALGORITHMS FOR CDMA SYSTEMS USING 
PILOT STRENGTH MEASUREMENTS, the specification of which is attached hereto. 

I hereby state that I have reviewed and understand the contents of the above identified 
specification, including the claims, as amended by an amendment, if any, specifically referred to 
in this oath or declaration. 

I acknowledge the duty to disclose all information known to me which is material to 
patentability as defined in Tide 37, Code of Federal Regulations, 1.56. 

I hereby claim foreign priority benefits under Title 35, United States Code, 119 of any 
foreign application^) for patent or inventor's certificate listed below and have also identified 
below any foreign application for patent or inventor's certificate having a filing date before that 
of the application on which priority is claimed: 

None 

I hereby claim the benefit under Tide 35, United States Code, 120 of any United States 
application (s) listed below and, insofar as the subject matter of each of the claims of this 
application is not disclosed in the prior United States application in the manner provided by 
the first paragraph of Title 35, United States Code, 112, I acknowledge the duty to disclose all 
information known to me to be material to patentability as defined in Title 37, Code of Federal 
Regulations, 1.56 which became available between the filing date of the prior application and 
the national or PCT international filing date of this application: 

None 

I hereby declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true; and further that these 
statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United 
States Code and that such willful false statements may jeopardize the validity of the application 
or any patent issued thereon. 
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receive the patent, and to transact all business in the Patent and Trademark Office connected 
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Lester H. Birnbaum (Reg. No. 25830) 

Richard J. Botos (Reg. No. 32016) 

Jeffery J. Brosemer (Reg. No. 36096) 

Kenneth M. Brown (Reg. No. 37590) 

Donald P. Dinella (Reg. No. 39961) 

Guy Eriksen (Reg. No. P-41736) 

Martin I. Finston (Reg. No. 31613) 

James H. Fox (Reg. No. 29379) 

William S. Francos (Reg. No. 38456) 

Barry H. Freedman (Reg. No. 26166) 

Julio A. Garceran (Reg. No. 37138) 

Mony R. Ghose (Reg. No. 38159) 

Jimmy Goo (Reg. No. 36528) 

Anthony Grillo (Reg. No. 36535) 

Stephen M. Gurey (Reg. No. 27336) 

John M. Harman (Reg. No. 38173) 

Donald E. Hayes, Jr. (Reg. No. 33245) 

John W. Hayes (Reg. No. 33900) 

Michael B. Johannesen (Reg. No. 35557) 

Mark A. Kurisko (Reg. No. 38944) 

Irena Lager (Reg. No. 39260) 

Christopher N. Malvone (Reg. No. 34866) 

Scott W. McLellan (Reg. No. 30776) 

Martin G. Meder (Reg. No. 34674) 

Geraldine Monteleone (Reg. No. 40097) 

John C Moran (Reg. No. 30782) 

Michael A. Morra (Reg. No. 28975) 

Gregory J. Murgia (Reg. No. 41209) 

Claude R. Narcisse (Reg. No. 38979) 

Joseph J. Opalach (Reg. No. 36229) 

Neil R. Ormos (Reg. No. 35309) 

Eugen E. Pacher (Reg. No. 29964) 
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Eli Weiss (Reg. No. 17765) 



I hereby appoint the attorney(s) on ATTACHMENT A as associate attorney(s) in the 
aforementioned application, with full power solely to prosecute said application, to make 
alterations and amendments therein, to receive the patent, and to transact all business in the 
Patent and Trademark Office connected with the prosecution of said application. No other 
powers are granted to such associate attorney(s) and such associate attorney(s) are specifically 
denied any power of substitution or revocation. 
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